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Abstract. —Advances in molecular techniques have facilitated use of genetic data in demographic wild¬ 
life studies. An important first step in genetic mark-recapture is selecting markers that uniquely “mark” 
and reliably “recapture” individuals. Markers should be tested on reliable DNAfrom known individuals 
(blood) before being used on non-invasively sampled DNA (hair, scat, or molted feathers). To evaluate 
whether Northern Goshawks (Accipiter gentilis ) can be uniquely identified by genotyping, 113 known 
(banded, sexed) goshawks from the Kaibab Plateau, Arizona, were genotyped using DNA from blood 
and five microsatellite markers and a sex-linked gene. We used mean relatedness to test whether adults 
in the population were related and probability of identity (P (ID) = probability that two random individ¬ 
uals from the population have the same genotype) to test the ability of multi-locus genotyping for 
uniquely identifying goshawks. We used genetic data to assess inbreeding and demographic data to 
estimate the effective population size. Sixty-nine adult goshawks were sexed correctly and genotyped. 
Expected heterozygosity was high (H £ = 0.81), and relatedness among adults was low (r = —0.017). All 
individuals sampled (69 adults, 44 nestlings) had unique five-locus genotypes, the overall probability of 
identity was low (P(id> unbiased ~ 7.03 X 10 -7 ), and the observed P^ was <0.0001. Thus, Kaibab goshawks 
were uniquely “marked” by genotyping. Despite a small effective population size (N e — 37 individuals), 
goshawks on the Kaibab Plateau functioned as a large breeding population with no inbreeding (F IS = 
—0.001). We hypothesized that genetic diversity is maintained by gene flow via immigration of individuals 
from distant forests. 

KEY WORDS: Northern Goshawk, Accipiter gentilis; capture-recapture, genetic tagging, individual identification, 
molecular sexing, probability of identity. 


GENETICA POBLACIONAL Y GENOTIPIFICACION PARA ESTUDIOS DE MARCADO-REGAPTURA 
DE ACCIPITER GENTILIS EN KAIBAB PLATEAU, ARIZONA 

RESUMEN. —Los avances en las tecnicas moleculares han facilitado el uso de informacion genetica en 
estudios demograficos de fauna silvestre. Un primer paso importante en estudios geneticos de marcado 
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y recaptura es seleccionar marcadores que “marquen” inequivocaiiiente y que permitan “recapturar” 
confiablemente a los individuos. Los marcadores deben ser probados en ADN confiable de individuos 
conocidos (sangre) antes de ser usados en muestreos no invasivos de ADN (pelo, excremento, o plumas 
mudadas). Para evaluar si los individuos de la especie (Accipiter gentilis) pueden ser idendficados por 
genotipificacion, estudiamos 113 gavilanes conocidos (sexados y anillados) de la meseta Kaibab, Arizona, 
usando ADN de la sangre y cinco marcadores microsatelites y un gen ligado al sexo. Usamos el par- 
entesco promedio para probar si los adultos en la poblacion estaban relacionados y la probabilidad de 
identidad (P(id> = la probabilidad que dos individuos al azar de la poblacion tengan el mismo genotipo) 
para probar la habilidad de la genotipificacion mediante multiples loci para identificar individualmente 
a los gavilanes. Utilizamos datos geneticos para determinar el grado de endogamia e informacion de- 
mografica para estimar el tamano efectivo de la poblacion. Sesenta y nueve gavilanes adultos fueron 
correctamente sexados y genotipificados. La heterocigocidad esperada fue alta (H e = 0.81) y el par- 
entesco entre adultos fue bajo (r = —0.017). Todos los individuos analizados (69 adultos, 44 polluelos) 
tuvieron genotipos unicos con cinco loci, la probabilidad de identidad total fue baja (P(id> unbiased = 
7.03 X 10 -7 ), y la P(id) observada fue <0.0001. Asi, los gavilanes de Kaibab fueron “marcados” singu- 
larmente por genotipificacion. A pesar de un tamano efectivo de la poblacion pequeno (N e = 37 
individuos), los gavilanes en el Kaibab Plateau funcionaron como una poblacion reproductiva grande y 
sin endogamia (Fjg = —0.001). Nuestra hipotesis es que la diversidad genetica es mantenida por flujo 
genetico a traves de la inmigracion de individuos de los bosques lejanos. 

[Traduccion de Mauricio Cotera] 


Northern Goshawks (Accipiter gentilis) are highly 
secretive and are most easily detected during the 
breeding season when they aggressively defend 
their nests and young. While their defensive be¬ 
havior at nests facilitates capture-recapture studies 
of breeding individuals, population monitoring is 
difficult because individuals often forgo breeding, 
or their nests fail early in a breeding season. Even 
in years of high productivity, mark-recapture stud¬ 
ies can be prohibitively expensive because popu¬ 
lation sampling requires large field crews and mul¬ 
tiple nest visits to many breeding territories to 
capture and recapture breeding goshawks (Reyn¬ 
olds et al. 2005). 

Because of recent improvements in molecular 
techniques (Haig 1998, Parkeuet al. 1998), genetic 
capture-recapture may be a viable alternative to tra¬ 
ditional capture-recapture methods for goshawks. 
Collecting molted feathers requires fewer nest vis¬ 
its than traditional capture-recapture methods. 
Breeding goshawks begin an annual molt during 
spring (Squires and Reynolds 1997) and because 
they spend much of the breeding season near their 
nests, they drop many of their molts within their 
nest areas, including years when nesting attempts 
fail. Thus, goshawk feathers are readily collected 
from nest areas and may provide an efficient 
means to non-invasively sample their populations. 

Several factors influence the success of genetic 
capture-recapture studies. An appropriate number 
of highly variable genetic markers for identifying 
individuals are required, and potential biases must 


be identified. “Shadow effects” (lack of discrimi¬ 
nation of individuals because of low variability or 
sampling too few markers) can negatively bias es¬ 
timates of population abundance and positively 
bias estimates of survival (Mills et al. 2000). On the 
other hand, when more markers than necessary 
are used, population abundance may be overesti¬ 
mated and survival underestimated if genotyping 
errors add unique “genotypes,” and thus individ¬ 
uals, to population samples (Lukacs and Burnham 
2005). Both biases will inflate variance and lower 
precision of parameter estimates (Lukacs and 
Burnham 2005). 

Microsatellites are currently a preferred molec¬ 
ular marker for identifying individuals because 
they are easily interpreted (i.e., heterozygous ge¬ 
notypes are easily distinguished from homozygous 
genotypes), highly variable, bi-parentally inherited, 
and generally appear to be selectively neutral. Fur¬ 
ther, a large body of literature exists on microsat¬ 
ellite evolution (Jarne and Lagoda 1996, Goldstein 
and Pollock 1997, Estoup et al. 2002), which has 
facilitated the development of much statistical the¬ 
ory and analytical software (Hedrick 2005). How¬ 
ever, microsatellites are expensive and time con¬ 
suming to develop for each newly-studied species. 
Occasionally primers used to amplify microsatellite 
markers in one species can be used in related spe¬ 
cies (Ellegren 1992, Primmer et al. 1996, Galbusera 
et al. 2000, Martinez-Cruz et al. 2002). 

Prior to starting a non-invasive genetic study, es¬ 
tablishing intrapopulation genetic structure (i.e., 
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levels of inbreeding and relatedness) and the fre¬ 
quency of null alleles (alleles that fail to amplify) 
is necessary for providing baselines against which 
feather samples can be compared (Mills et al. 
2000). Likewise, it is important to establish statis¬ 
tical power of multi-locus genotyping for identify¬ 
ing individuals with an independent population 
sample. We present results from a pilot study where 
we assessed the feasibility of implementing a non- 
invasive genetic capture-recapture study on a pop¬ 
ulation of Northern Goshawks on the Raibab Pla¬ 
teau, Arizona. Before assessing the utility of molted 
feathers as a viable source of DNA, we established 
a dependable genotyping marker set using DNA 
derived from blood (Taberlet and Luikart 1999). 

Our objectives were to: (1) screen species-specif¬ 
ic and cross-specific (among species) microsatellite 
markers, (2) test a sex-linked gene in goshawks 
known to distinguish males and females in other 
raptors (Kahn et al. 1998), (3) assess the power of 
multi-locus genotyping to uniquely identify individ¬ 
uals using probability of identity analysis (P (ID) ; 
probability that two individuals drawn at random 
from the same population share the same multi¬ 
locus genotype), and (4) estimate average relat¬ 
edness, inbreeding, and effective population size 
for the goshawk population on the Kaibab Plateau. 

Methods 

Field Collection. The goshawk study population is lo¬ 
cated on the Kaibab Plateau in northern Arizona, an area 
that includes the North Kaibab Ranger District of the 
Kaibab National Forest and the North Rim of the Grand 
Canyon National Park (for descriptions of the study area 
see Reich et al. 2004, Reynolds and Joy 2005). It is a 
forested plateau surrounded by shrub-steppe habitat— 
the nearest forests being 97 km to the north, 250 km to 
the east, 80 km to the west, and 89 km to the south— 
except for a small patch of forest 18 km south on the 
south rim of the Grand Canyon. Sampled nests were well 
distributed across the study area. We captured 69 adult, 
breeding goshawks (1991-93, 2000-02) and 44 of their 
nestlings (Reynolds et al. 1994). We sexed adult goshawks 
using morphometries (mass, tarsus length) and behavior. 
Blood was sampled from the brachial vein with 22-gauge 
needles and non-heparinized capillary tubes (volume ^ 
0.10 ml). Blood was transferred into STE (Sodium Chlo- 
nde-Tris-EDTA) buffer-filled storage tubes kept cool in 
insulated containers with frozen cold-packs until crews 
returned to the field station, where samples were subse¬ 
quently frozen (—20°C). At the close of the field season, 
blood was transferred to and stored at — 80°C at Colorado 
State University, Fort Collins, CO. 

Laboratory Methods. We extracted DNA using 
QIAamp mini blood kits (Qiagen, Inc., Valencia, CA 
U.S.A.) following the manufacture’s protocol. To find mi¬ 
crosatellites, we screened published and unpublished 


primer sets that included microsatellites originally isolat¬ 
ed from Northern Goshawks (Topinka and May 2004), 
European goshawks ( Accipiter gentilis gentilis ), Golden Ea¬ 
gles ( Aquila chrysaetos), and Red Kites ( Milvus milvus, 
Peck 2000). We also tested primers that amplify an intron 
within the avian CHD (chromo-helicase-DNA binding), 
which was used to determine gender in Red-tailed Hawks 
(Buteo jamaicensis) and Great Horned Owls ( Bubo virgim- 
anusr, Kahn et al. 1998). The CHD gene is located on the 
Z and W sex chromosomes. We expected males to be 
homozygous (ZZ genotype) and females to be heterozy¬ 
gous (ZW genotype). 

We used PCR (Polymerase-Chain-Reaction) to amplify 
microsatellites in 25 |xl reactions using 0.5 pi (AGE la) 
or 1.0 pi (all other markers) of template DNA, 2.5 pi 
10X buffer containing 15 raM Mg 2 Cl 3 (Promega Corp. 
Madison, WI U.S.A.; for markers AGE la, AGE 2 and 
AGE 4 an additional 3 mM Mg 2 Cl 3 was added), 20 mM 
dNTPs, 25 pM each primer, 1U Taq polymerase, and one 
drop mineral oil to prevent evaporation. Negative con¬ 
trols (reactions that include all reagents except template 
DNA) were included in every set of reactions, and we 
used “cold start” PCR where tubes (in racks) were kept 
on ice to prevent premature non-specific priming. We 
used MJR PTC-100 thermocyclers programmed for the 
following protocol: denature at 94°C for 4 min, 31 cycles 
of denature at 94°C for 40 sec, annealing at 58°C for 40 
sec, and chain extension at 72°C for 40 sec, with a final 
extension at 72°C for 5 min. 

We used PCR to amplify the CHD sex-linked gene in 
25pi reactions using l.Opl template DNA and the same 
reaction buffer described above. The PCR protocol in¬ 
cluded an initial 5 min at 95°C denature, 11 cycles of 
denature at 94°C for 30 sec, annealing at 52°C for 35 sec, 
and chain extension at 72°C for 2.0 min, 31 cycles of 
denature at 92°C for 30 sec, annealing at 56°C for 35 sec, 
and chain extension at 72°C for 2.0 min, with a final 
extension at 72°C for 7 min. 

We used gel electrophoresis to separate alleles. For mi- 
crosatellites, we used 8% polyacrylamide (Long Ranger, 
Cambrex Corp., Rockland, MA U.S.A.) denaturing gels 
(55 cm long) that were run at 45 watts for 4-5 hr, de¬ 
pending on allele size. For the CHD sex-linked gene, we 
used single-strand-conformation-polymorphism (SSCP) 
methods (Hiss et al. 1994) and electrophoresed alleles 
on non-denaturing gels at 5 watts for 15 hr. 

For microsatellites, we established allele standards us¬ 
ing representative samples from our first gel and then 
standardized all other gels using those same samples. 
Gels were scored visually, and allele standards were run 
on both sides of a gel to account for gel ambiguities that 
cause slight variations in migration distances. Further, a 
subset of individuals (N = 23) was genotyped a second 
time to validate scores for microsatellite markers. For the 
CHD marker, we ran all known females together (N = 
40) and all known males (N = 29) together to familiarize 
ourselves with allele morphology. Although not se¬ 
quenced, the fragments were ca. 240-260 base pairs, and 
Z and W alleles were similar in size, but were differenti¬ 
ated by the SSCP analysis (Hiss et al. 1994). 

Population Genetic Analysis. Population substructure, 
inbreeding, and genetic drift can reduce heterozygosity 
in populations. However, low yield and degraded DNA 
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sampled from sources such as molted feathers can arti¬ 
ficially reduce population heterozygosity if allelic drop¬ 
out (ADO; one of two alleles in a heterozygous individual 
fails to amplify) at one or more markers occurs. It is, 
therefore, important to use high yield sources of DNA 
(typically blood) from a known reference population to 
determine frequency of ADO (or null alleles) and true 
levels of heterozygosity (Taberlet et al. 1999). 

We used Cervus 2.0 (Marshall et al. 1998) to estimate 
observed (H 0 ) and expected (H E ) population heterozy¬ 
gosity and null allele frequencies. Cervus provides esti¬ 
mates of null allele frequencies with an iterative algo¬ 
rithm based on differences between observed and 
expected homozygote frequencies. We used Genepop 3.4 
(Raymond and Rousset 1995) to test for departures from 
Hardy-Weinberg equilibrium (random mating) and GDA 
1.0 (Lewis and Zaykin 2001) to test for linkage-disequi¬ 
librium (genotypes at one marker are independent from 
genotypes at other markers) and to estimate F IS , an in¬ 
dicator of population substructure and inbreeding. For a 
review of /-statistics and microsatellite genetic markers, 
see Balloux and Lugon-Moulin (2002). 

To test our assumption that our sample of adult gos¬ 
hawks was not comprised of closely-related individuals, 
we used Identix 1.1 (Belkhir et al. 2002) to estimate 
mean pairwise relatedness. We used Queller and Good¬ 
night’s (1989) estimator option, and tested the null hy¬ 
pothesis of no relatedness by comparing our estimate to 
a distribution of coefficients derived through convention¬ 
al Monte Carlo resampling procedures (1000 permuta¬ 
tions). 

Probability of Identity. The uniqueness of an individ¬ 
ual’s genotype depends on the number and polymor¬ 
phism (heterozygous) of the markers. Multi-locus geno¬ 
types based on few highly-variable markers can be as 
powerful as those based on many less variable markers 
(Waits et al. 2001). Mills et al. (2000) suggested for stud¬ 
ies of genetic demography that profiles should be based 
on multi-locus genotypes capable of discriminating indi¬ 
viduals with 99% certainty. Estimating probability of iden¬ 
tity (P(id)) is one way to establish this certainty when it is 
expressed as 1 — P (ID) . P(id) is similar to the match prob¬ 
ability used in human forensics (Evett and Weir 1998, 
Avise 2004, Hedrick 2005), but is less susceptible to vio¬ 
lations of linkage-disequilibrium and Hardy-Weinberg 
equilibrium, both of which can be prevalent in small, iso¬ 
lated, or substructured populations (Waits et al. 2001). 

P (1D ) analysis includes two steps. First, two theoretical 
P (ID )’s, one for unrelated individuals (P(id) unbiased)) an d 
one for siblings (P ( i D ) sibs ), are estimated (for equations 
see Waits et al. 2001). Both estimators use population 
allele frequency data, and P(iD)unbia S ed* s corrected for bias 
in small samples. The two estimators provide lower and 
upper confidence bounds on the number of markers 
needed to discriminate individuals accurately. If the study 
population is composed of many related individuals, then 
resolving those individuals requires more markers. Step 
two involves calculating an observed P ( iD)obs based on ac¬ 
tual multi-locus genotypes from a known population sam¬ 
ple and is simply the proportion of all possible pairs of 
individuals with identical multi-locus genotypes (Waits et 
al. 2001). 

To estimate both theoretical P (ID) ’s and to quantify 


P(iD)obs (the proportion of individuals that share geno¬ 
types), we used PROB-ID5 (Waits et al. 2001). We used 
multi-locus genotypes derived from 69 adult goshawks, 
which we assume to be unrelated (see below), and 44 of 
their nestlings (sibling groups of 2-4 nestlings). We first 
analyzed the adults and then added the offspring/sibling 
groups. We used all five microsatellite markers and the 
CHD sex-linked gene and added markers sequentially 
starting with those having the highest number of alleles. 

Effective Population Size. To evaluate whether immi¬ 
gration and gene flow influenced genetic structure of the 
Kaibab population we estimated its effective size (N e ) 
Effective population size is the idealized number of in¬ 
dividuals in a population measured either demographi- 
cally, in terms how many individuals actually contribute 
to breeding (i.e., variance in productivity), or genetically, 
using f^-statistics and measures of inbreeding where the 
assumption of non-overlapping generations exists (Bar¬ 
ton and Whitlock 1997, Hedrick 2005). Goshawks sam¬ 
pled on the Kaibab Plateau during the study period likely 
represent at least three overlapping generations; thus, we 
relied on reproductive data to estimate N e . 

Effective population size is generally smaller than the 
censused population ( N ). Counts of breeding pairs of 
goshawks can be used to index N e , but not all goshawk 
pairs on the Kaibab produced an equal number of off¬ 
spring during the study (Wiens and Reynolds 2005). We 
therefore estimated annual N e ’s (equation 6.8a in Hed¬ 
rick 2005) as: 


where k is mean productivity measured as the number of 
young fledged per used nest, eggs laid, subsequent young 
fledged, or the nest failed (1991-2003; Reynolds et al. 
2005), V k is the variance in annual mean productivity, and 
N is the annual count of breeding pairs for the year. We 
then calculated a 13-yr harmonic mean of annual iVe’s 
(equation 6.12b in Hedrick 2005) for our final size esti¬ 
mate. 

Results 

Genetic Markers. Of nine cross-specific and sub¬ 
specific markers tested, two did not resolve alleles 
successfully, six amplified successfully but lacked 
variability, and one both amplified and was poly¬ 
morphic (AGE la, Table 1). All four microsatellites 
originally isolated from Northern Goshawks ampli¬ 
fied and were polymorphic (Table 1). The CHD 
sex-linked gene amplified and SSCP genotypes 
were consistent within the sexes (females, N = 40; 
males N = 29), making it useful for distinguishing 
between male and female goshawks. We validated 
our amplification and scoring of microsatellite 
markers after all individuals were genotyped and 
scored the first time. We genotyped the 23 individ¬ 
uals used as standards a second time using DNA 
that was archived and remained untouched in our 



Table 1. Statistics for microsatellites tested on blood-derived DNA sampled from female (N = 40) and male (N — 29) Northern Goshawks ( Acapitergentilis), 
Kaibab Plateau, AZ (1991-1993, 2000-2002). 
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freezers from the time of original DNA extraction. 
We found only a single scoring error out of 230 
opportunities (23 samples genotyped twice for five 
markers). This was a recording error; the sample 
had actually genotyped correctly. 

In comparisons of expected (H £ ) and observed 
(H 0 ) heterozygosity at each marker, four of the five 
microsatellites were similar. However, one marker 
(AGE 6) significantly departed from Hardy-Wein- 
berg expectations (P < 0.01) due to a deficiency 
of heterozygote genotypes (Table 1). Based on our 
data, we suspected this marker was sex-linked, as 
we found strong linkage-disequilibrium between 
AGE 6 and the sex-linked CHD marker (P< 0.001, 
Fisher’s method, 3200 runs) indicating that the 
two markers segregate together. It appeared that 
the marker was located on the Z sex-chromosome 
because all females (ZW) had only a single allele 
(homozygous), while most males (ZZ) were hetero¬ 
zygous. We hypothesize that alleles on the female’s 
W-chromosome are non-amplifiable (null) because 
of mutations in the priming sequences flanking the 
marker (Scribner and Pearce 2000) or because the 
marker on the Z-chromosome simply has no ho¬ 
mologous region on the W-chromosome. We 
found no evidence for linkage-disequilibrium be¬ 
tween the other four microsatellite markers (Table 
1) when AGE 6 and CHD were excluded from the 
analysis. We found no evidence of null alleles, 
which is important for future assessments of ge- 
notyping error when using feathers as a source of 
DNA. 

Population Genetics. F IS measures departures of 
observed and expected heterozygosity under as¬ 
sumptions of random mating and indicates either 
inbreeding (F IS > zero) or inbreeding avoidance 
(F IS < zero). Thus, highly structured or isolated 
populations that experience genetic drift generally 
exhibit positive F IS values. Alternatively, large pop¬ 
ulations or those experiencing high gene flow gen¬ 
erally exhibit nonsignificant or negative F IS values. 

We found no evidence of inbreeding or inbreed¬ 
ing avoidance (Fj S = —0.001; 95% Cl = —0.070- 
0.063; AGE 6 excluded; Table 1), suggesting that 
Kaibab goshawks mate randomly. Lack of inbreed¬ 
ing could result from large population size, gene 
flow by immigrants, or both. However, our demo- 
graphically derived estimate of effective population 
size (N e = 37 individuals; range = 10-86) indicated 
that the population was demographically small, 
thus making gene flow a more likely source of ge¬ 
netic variability. This is consistent with our estimate 
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Mean Coefficient of Relatedness 

Figure 1. Observed mean relatedness relative to 1000 
randomized populations assumed to lack relatedness. 
The observed mean falls below that which is expected at 
random, occurring with a probability of 3.7%, indicating 
that Northern Goshawks (Accipiter gentilis ) on the Kaibab 
Plateau, Arizona, are less related than expected at ran¬ 
dom. 


of relatedness among adults goshawks, where mean 
relatedness ( r^ = —0.017) was less than expected 
by random (Fig. 1). 

Probability of Identity. Certainty of individual 
identification is equal to 1-P (ID) , and therefore the 
goal in estimating probability of identity is to ob¬ 
tain small values of P (ID) . Waits et al. (2001) sug¬ 
gested a value ^0.0001 for forensic investigations 
where estimates of demographic parameters are 
needed. This threshold is interpreted as a 1:10 000 
chance that two individuals sampled from the same 
population will have the same multi-locus geno¬ 
type. 

We found that all 69 unrelated adults had 
unique multi-locus genotypes with the inclusion of 
the first three markers (P ( i D)0 b s < 0.0001), and like¬ 
wise the estimated P(n» met the 0.0001 threshold 
(P( ID) unbiased = 1*13 X 10~ 4 ; Fig. 2a). With five 
markers, the same sample had a P(iD) U nbiased = 7.03 
X 10 -7 . However, based on demographic data, we 
know that siblings and parent-offspring nested si¬ 
multaneously in the Kaibab population (R. Reyn¬ 
olds unpubl. data). To model this effect we added 
44 nestlings-siblings to the sample. While the two 
theoretical P( ID) ’s did not change, all five markers 
were required to differentiate individuals (P(iD)obs 
< 0.0001, Fig. 2b). 

In both cases (adult only and adults with off¬ 
spring-sibling groups), our sample of markers was 
insufficiently large to bring the P(iD) s ibs to the 
0.0001 threshold. Thus, we were not able to esti¬ 
mate an upper number of markers needed for this 
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Figure 2. Relationship between theoretical, observed, 
and sib probability of identity (P (ID) ) for Northern Gos¬ 
hawks (Accipiter gentilis) on the Kaibab Plateau, Arizona. 
The first five markers are microsatellites and the sixth is 
the CHD sex-linked gene. Observed data closely tracked 
that of the theoretical estimator; however, while (a) all 
69 unrelated adult goshawks were resolved after the first 
three markers (observed P (ID) < 0.0001), (b) it took an 
additional two markers to resolve sibling and parent-off¬ 
spring goshawks when 44 nestlings were added to the 
sample. In both cases, the theoretical P (ID) met our 
0.0001 threshold (a 1:10 000 chance that two individuals 
sampled from the same population have identical multi¬ 
locus genotypes). However, we did not analyze enough 
markers to bring the sibling P(idj to the threshold level. 


resolution. Nonetheless, with five markers P ( iD) S ibs 
= 6.17 X 10 3 , which translates into a six in a 1000 
chance of drawing two identical genotypes. Be¬ 
cause we sampled many parent-offspring pairs that 
we could nevertheless distinguish, we are confident 
the five combined markers provide unique genetic 
marks. 

Discussion 

Our intent in this study was to develop a set of 
genetic markers that uniquely identified individual 
goshawks. We desired to establish this marker set 
using high-yield DNA (blood) sampled from a 
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known reference population. While most microsat¬ 
ellites tested did not amplify or were monomor- 
phic (most A. g. gentilis markers), we did find a set 
of highly variable markers that consistently ampli¬ 
fied DNA from blood. Elsewhere (S. Bayard de 
Volo unpubl. data), we found that the same genetic 
markers consistently and reliably amplified DNA 
from molted feathers. We note that because AGE 
6 is probably on the Z sex chromosome, its utility 
for estimating within population relatedness and 
levels of inbreeding is limited to samples from 
males. However, differences in allele frequencies 
between populations will still be useful for larger- 
scale studies comparing populations. In a study of 
goshawks in Utah, Sonsthagen (2002) used a dif¬ 
ferent set of microsatellite markers than ours. Of 
eight markers, only one of theirs exhibited the 
same number of alleles (11 alleles, H 0 = 0.73, H £ 
= 0.74) as AGE 6 did in our study. This alternative 
marker (BV 20; Gautschi et al. 2000) would be use¬ 
ful if it exhibited similar levels of heterozygosity 
and allelic diversity in the Kaibab population. Re¬ 
placing AGE 6 with a less variable marker would 
result in having to add more markers to the entire 
genetic profile, which would introduce more op¬ 
portunities for genotyping error. We are currently 
testing BV 20 to see if it is an effective replacement 
for AGE 6. 

With the five microsatellite markers tested (Ta¬ 
ble 1), all 113 goshawks sampled had unique multi¬ 
locus genotypes resulting in a P (ID)obs < 0.0001 
(Fig. 2b) and a P ( i D ) U nbiased = 7.03 X lO -7 . This was 
a powerful result considering that our sample in¬ 
cluded many parent-offspring and sibling pairs 
from the same nest. Likewise, the five microsatel¬ 
lites showed a high level of expected heterozygosity 
(H e = 0.81). Others have shown that marker sets 
composed of five markers that result in H E s 0.80 
will have a theoretical P (ID) ^ 0.0001 (Waits et al. 
2001). In Paetkau’s (2003) retrospective analysis of 
21 non-invasive genetic studies in bears (Ursus 
americanus, U. arctos ), the number of markers used 
was determined by whether the first five most-var¬ 
iable microsatellite markers together had H E > 
0.80. He found that for some black bear popula¬ 
tions H e was >0.80 for five markers; however, for 
others, and for all grizzly bear populations, that H E 
was <0.80, requiring the marker set to be in¬ 
creased to six or seven loci in order to discriminate 
among individuals. 

We note that our estimates of P< ID) are specific 
to the Kaibab goshawk population; we cannot pre¬ 


dict with complete certainty that these same mark¬ 
ers will uniquely mark goshawks from other pop¬ 
ulations. Power of discrimination depends on 
population-specific levels of genetic variability (het¬ 
erozygosity) ; goshawk populations that are less var¬ 
iable because of geographic isolation or habitat 
fragmentation may require more markers to 
uniquely genotype individuals (Paetkau and Stro- 
beck 1994). However, goshawks are highly vagile, 
and we suspect gene flow is high among popula¬ 
tions. These goshawk populations will probably ex¬ 
hibit similar heterozygosity, and the marker set test¬ 
ed here should prove useful for other studies. 

The Kaibab goshawk population exhibits high 
genetic variability (Table 1), despite its geographic 
isolation and small effective population size (based 
on demographic data; 13-yr x = 37 individuals). 
Several explanations may account for this. First, it 
is possible that the markers used in this study are 
under selective sweeps with genes that are affected 
by balancing selection for heterozygous genotypes. 
Such selection has been found for the genes of the 
MHC (major histocompatability complex) in mam¬ 
mals (Avise 2004), in which heterozygous individ¬ 
uals experience a fitness advantage. However, we 
suggest it is unlikely that all four nonsex-linked 
markers would be under the same selective pres¬ 
sures, given that they exhibit independent segre¬ 
gation (no evidence of linkage disequilibrium; see 
Black et al. 2001). 

A second and more likely explanation is that ac¬ 
tual N e for this goshawk population is much larger 
because geographically distant populations in the 
region are connected by migration and gene flow. 
While adult goshawks are mostly sedentary on 
breeding territories (Detrich and Woodbridge 
1994, Squires and Ruggiero 1995, Reynolds and Joy 
2005), band recoveries of first-year goshawks from 
the Kaibab Plateau indicate dispersal distances of 
up to 440 km (Wiens 2004). In addition, telemetry 
data show that juvenile goshawks disperse from the 
Kaibab Plateau in their first year, with the majority 
moving beyond the 80 km detection distance 
(Wiens 2004). Further, Wiens (2004) showed that 
only 11% of 614 banded nestlings returned to be 
recruited into the Kaibab breeding population, in¬ 
dicating high first-year mortality or low natal-site 
fidelity. Evidence for the latter is indicated by the 
lack of population genetic structure for goshawks 
in Utah (Sonsthagen et al. 2004). To better assess 
actual effective population size for goshawks in 
western North America, we are expanding our 
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studies to include populations in the western por¬ 
tion of the species range. Data from these studies 
should allow a more comprehensive evaluation of 
the genetic structure and effective population size 
for goshawks in the West. 

Conclusions 

Genetic marking of Northern Goshawks on the 
Kaibab Plateau is both feasible and reliable. Like¬ 
wise, non-invasive genetic sampling will provide an 
alternative method for demographic and genetic 
data collection, as we have found that molted 
feathers are as reliable a source of DNA as blood 
(S. Bayard de Volo unpubl. data). Because gos¬ 
hawks show high territory fidelity (e.g., Detrich 
and Woodbridge 1994, Reynolds and Joy 2005), 
they are particularly well suited for non-invasive ge¬ 
netic sampling. We recommend that monitoring 
programs implement rigorous field collection of 
molted feathers. As with any demographic study, 
valid inferences to the population depend on ap¬ 
propriate spatial and temporal sampling from that 
population. Researchers and managers interested 
in implementing non-invasive genetic mark-recap¬ 
ture to study goshawks should contact the corre¬ 
sponding author or refer to Bayard de Volo (2005). 
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Appendix. Source, repeat qualities and accession numbers or primer sequences for micro satellites found to be useful 
for genotyping Northern Goshawks {Accipiter gentilis) on the Kaibab Plateau, Arizona. 


Marker 

Repeat Motif 3 

Allele Size 
in Base 

Pairs 

Accession No. 
or Primer Sequences 

5' to 3' 

Author 

AGE 1 

(gggaa) 8 . . (gaga) 9 . . (gagaa) 3 

216 

AY312451 


AGE 2 

(gagaa) 10 (ga) 4 

170 

AY312452 

Topinka and 

AGE 6 

(gagaa) 4 . .(gagaa) 2 . . (gagaa) 5 

259 

AY312456 

May 2004 

AGE 4 

(gagaa) 19 

275 

AY312454 


AGE la 

(ggat) s 

208 b 

f acaactgggctgtgctttgc 
r cttcccggtggctgaggctt 

Peck 2000 


a Sequenced by authors. 

b Mean allele size in European goshawk (A. g. gentilis). 





